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Abstract 

A novel expansion of the evolution operator associated with a - in general, time-dependent 
- perturbed quantum Hamiltonian is presented. It is shown that it has a wide range of 
possible realizations that can be fitted according to computational convenience or to satisfy 
specific requirements. As a remarkable example, the quantum Hamiltonian describing a 
laser-driven trapped ion is studied in detail. 
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1 Introduction 



The explicit determination of the evolution operator associated with a quantum system — 
namely, the determination of its explicit action on the state vectors, or, equivalently, on 
a given orthonormal basis — is, in general, a 'touchy business'. If the Hamiltonian of the 
system does not depend on time and has the form of the sum of a solvable unperturbed 
Hamiltonian plus an analytic perturbation, one can use the tools of standard perturbation 
theory for linear operators, based on the expansion of the resolvent, in order to get 
approximate expressions of the evolution operator. Alternatively, one can apply a suit- 
able operator perturbative approach [21 El EJ El E] in order to obtain a very convenient 
perturbative expansion of the evolution operator whose truncations have the remarkable 
property of forming one-parameter groups of unitary operators. The main goal of this 
paper is to extend, in a natural way, such a perturbative approach to the general case 
where the quantum Hamiltonian may depend on time. 

* Paolo.Aniello@na.infn.it 
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We stress that if a quantum Hamiltonian is time-dependent — i.e. it describes a non- 
isolated quantum system — the task of determining the associated time evolution operator 
is a tough problem since, whenever the values of the Hamiltonian at different times do not 
commute, the evolution operator does not admit a simple formal expression. 
In two fundamental papers 0, Dyson developed an expansion of the evolution operator 
that has been adopted extensively in any field of physics. Dyson expansion has a trans- 
parent physical interpretation in terms of time ordered elementary processes which makes 
its application particularly appealing, especially in quantum field theory. On the other 
hand, for many applications, Dyson expansion has severe drawbacks, as a low convergence 
rate and the lack of unitarity of its truncations [S] . 

Later, Magnus [Oj introduced an expansion of the evolution operator such that each of 
its truncations retains the property of being unitary. Magnus expansion has been 're- 
discovered' and re-elaborated several times, and applied successfully to several problems 
(see |10j and references therein). It is written in the form of the exponential of the expan- 
sion of a suitable time-dependent anti-hermitian operator which can be deduced, order 
by order, from the Hamiltonian of the system. Now, precisely for this reason — just like 
for the evolution operator generated by a time-independent Hamiltonian — the problem 
of computing explicitly the action of (any truncation of) the Magnus expansion on the 
state vectors is non-trivial. Truncating the power expansion of the exponential would 
lead to non-unitary expressions, thus to the loss of the most important feature of Magnus 
expansion. Besides, in the important case of a time-independent Hamiltonian, there is 
no link between Rayleigh-Schrodinger-Kato perturbation theory for linear operators and 
Magnus expansion, a clue suggesting that the power of the perturbative approach is not 
'fully exploited' by this expansion. 

Then, the issue of finding a generalization of Magnus expansion retaining the property of 
having unitary truncations, but allowing more convenient solutions — or solutions hav- 
ing special properties — arises in a natural way. In the present paper, we have tried to 
achieve this result. We introduce a perturbative decomposition of the evolution operator 
that generalizes Magnus expansion, and opens the possibility of obtaining computational 
advantages and of satisfying specific requirements for the perturbative solutions by suit- 
ably fixing, at each perturbative order, certain arbitrary operators or operator-valued 
functions. It is important to remark that, in the case where the Hamiltonian does not 
depend on time, our decomposition, differently from Magnus expansion, assumes a special 
meaning and has a precise link with standard perturbation theory for linear operators. 
The idea of generalizing Magnus expansion is not completely new. It appears in a paper 
by Casas et al. in which the authors introduce the Floquet- Magnus expansion for the 
evolution operator associated with a (interaction picture) Hamiltonian depending period- 
ically on time. However, it turns out that our approach generalizes the one proposed by 
Casas et al. even in the case where the interaction picture Hamiltonian is periodic on time. 

We have made the choice of skipping, as far as possible, mathematical complications. 
Our choice is motivated by various reasons. First, we believe that heuristic investigation 
should always precede rigorous re-elaboration. Once that it is clear what the basic 'rules 
of the game' are, one can adopt the most appropriate mathematical tools. Moreover, we 
avoid the risk of hiding in a thick cloud of technicalities the main ideas and of discouraging 
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those who may want to apply our method for solving problems. It should be also observed 
that a recent trend in quantum mechanics is to focus on systems which can be described 
by effective Hamiltonians in finite-dimensional Hilbert spaces (consider the huge research 
area related to quantum computation and quantum information theory; see ref. ^2] an d 
the rich bibliography therein). The study of these systems is not affected by the techni- 
calities associated with the infinite-dimensional spaces but retains all the most intriguing 
features of quantum physics. 

The structure of the paper is the following. We begin introducing, in sect. [2 before going 
through the details of a our perturbative approach, a significant example — the quantum 
Hamiltonian describing a laser-driven trapped ion — example that will play a double role: 
a starting motivation for our analysis and a concrete playing field for practising our theo- 
retical results. In sect. |3] we establish the general form of the perturbative decomposition 
of the evolution operator that we propose. In the subsequent two sections, we pursue the 
task of finding a recursive procedure that allows to compute, order by order, the various 
terms of our perturbative expansion. Precisely: in sect. 12] we study the case where the 
Hamiltonian does not depend on time and investigate the link with standard perturbation 
theory for linear operators; in sect. [5] we consider the general case of a time-dependent 
Hamiltonian. Next, sect. El is devoted to study in detail the significant example introduced 
in sect. in such a way to illustrate the main features of our perturbative expansion of 
the evolution operator. Finally, in sect. conclusions are drawn, with a quick glance to 
further applications. 

2 A remarkable example: the ion trap Hamiltonian 

Before introducing our perturbative approach, we think that is worth considering first a 
remarkable example: the Hamiltonian describing a one-dimensional laser-driven ion trap. 
The study of this Hamiltonian, which will be performed systematically in sect. |H1 allows 
to illustrate in a simple way all the main points of the theory developed in the subsequent 
sections. 

A two-level ion of mass /x in a potential trap, with strong confinement along the y 
and z axes, and weak harmonic binding of frequency v along the x-axis (the 'trap axis'), 
can be described — neglecting the motion of the ions transverse to the trap axis — by a 
Hamiltonian of the following type (K= 1): 



where n = a^a is the 'number operator' — with a denoting the vibrational annihilation 
operator 



- and a z the effective spin operator associated with the internal degrees of freedom of 
the ion. Let us suppose now that the ion is addressed by a laser beam of frequency a in 



Hq = v h + — ea. 
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a 'traveling wave configuration'. Then, the Hamiltonian describing the physical system 
becomes (see, for instance, ref. |T3]): 

Ho + H t (t), (1) 
where the time-dependent interaction term H^(t) is defined by 

H t (t) := n R (e iat D(irtf a_ + e~ iat D(irj) a + ) , (2) 



with Q R denoting the Rabi frequency (which is proportional to the intensity of the laser 
field) and with a± := |±) (=p| , cr+ = o\L (a z |±) = ± |±)). Moreover, we have set: 



D(ii]) := exp \ irj ya + a ^jj j (3) 

where 

^ _ fcx cos ^ 

\J1\iV 

— with the wavevector and eft the angle between the x-axis and Icl — is the so-called 
'Lamb-Dicke parameter'. In the case where ?) <C 1 ('Lamb-Dicke regime') — a case often 
occurring in applications 1 — one can keep only those terms in the power expansion of 
D(irj) which are at most linear in r/: 

H l {t)^VL R (e iat {l-ir 1 {a + cJ))<T„+ h.c}j =: H$(t) . (5) 

Observe that the problem of dealing with the time-dependent Hamiltonian Hq + H^(t) 
can be bypassed by switching to the interaction picture with reference Hamiltonian ^aa z . 
Indeed, setting 

Rt := exp (-^aa z t\, (6) 
one obtains the time-independent 'rotating frame Hamiltonian' 

R\ (h + H : (t)-^aa z ^j R t = uh + ^5a z + n R [d^Y o_ + D(irj) 04 
(Lamb-Dicke regime: rj <C 1) ~ £Ir ( (1 — ir](a + a')) cr_ + h.c. ) , (7) 



where 5 := e — a is the ion-laser detuning. In many applications, the condition £l R <C is 
also verified; hence, it is natural to introduce the (dimensionless) perturbative parameter 

A^«l. (8) 



1 Notice that in a trap with a linear geometry, like the one we are considering here, one can modify the 
incidence angle <f) of the laser beam with respect to the principal axis of the trap in order to control the 
Lamb-Dicke parameter n. 
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Thus, the study of the time-dependent Hamiltonian can ultimately be reduced to the 
study of the simpler time-independent Hamiltonian 

H :=un + ^Sa z + Az; ((1 - 477(0 + a*)) a _ + h.c)j , (9) 

which has the form of a trivially solvable Hamiltonian 2 Hq := vn + \ ba z plus a 'small 
perturbation' := Xv ((1 — irj(a + a')) <r_ + /i.e.) . 

The exact eigenvalues and eigenprojectors of the Hamiltonian H, despite its formal 
simplicity, are not known. Hence, at this point, it would seem reasonable to adopt a 
(time- independent) perturbative approach for studying the Hamiltonian H. However, 
this is not what it is usually done in the literature. In fact, with the aim of providing an 
approximate expression of the evolution operator associated with the Hamiltonian H, the 
rotating wave approximation (RWA) — see, for instance, ref. ^1] — is usually applied. 
This approximation amounts to passing to a further interaction picture with reference 
Hamiltonian Hq, so obtaining the new interaction picture Hamiltonian 

H int (t) =\v(l- irj(e iut a t + e~ ivt a)) e~ iSt a_ + h.c. 

- which is once again time-dependent like the ion trap Hamiltonian ((J) — and, then, 
retaining only those terms in Hi nt (t) which are 'slowly rotating'; all the other terms, often 
called 'counter-rotating terms' (for historical reasons), are simply ignored. 
In particular, in correspondence to the three types of resonance condition 

5 = e-a&0, 5 + S-umO, (10) 

one obtains, respectively, the following three types of effective interaction picture Hamil- 
tonian: 

(f«Q) (11) 

(*«-!/) (12) 
(*«+!/) (13) 

These effective Hamiltonians, in correspondence to the respective resonances, commute 
with the reference Hamiltonian Hq. This is due to the fact that the resonances (|10|) are 
associated with the appearance of degeneracies in the spectrum of the reference Hamil- 
tonian, and the degenerate eigenspaces of Hq are invariant subspaces for the effective 
Hamiltonians Hz*, HQ, H^l, respectively 3 (in fact, it turns out that the spectrum of 
Ho is degenerate if and only if the condition |<5| = mis, m = 0,1,2, ... , holds). As a 
consequence, the evolution operators generated by the effective Hamiltonians H^, H^, 

2 Namely, a Hamiltonian such that a complete orthonormal set of eigenvectors is given simply by the 
standard basis {|n} <g> |±)}. 

3 In correspondence to the resonances 11UI . the only non-degenerate eigenspaces of Ho are given by 
span(jO) ® |+)), for 6 = — u, and span(|0) ® |— )), for 5 = v, that are invariant subspaces for and ffjg>, 
respectively. All the other eigenspaces are doubly degenerate. 



n eft 

fr(-) 
^eff 

n eff 



\v (cr_ + (T+) , 
iXvrj (a) a + — a a_ ) , 

i \ v i] \aa + — a^ O- 
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i/gg can be explicitly determined. 

Notice moreover that, in particular, H^l — up to a unitary transformation 

gig g-iTrn/2 = A 1/ 7/ (a <7+ + 0* <7_) 

— has the same form of the typical interaction term of the Jaynes-Cummings Hamiltonian 
(see ref. [TK]). 

At this point, it is quite natural to address the following questions: 

1. Is it really necessary, after having obtained the time-independent Hamiltonian H, 
to switch to a further interaction picture in order to achieve a 'good approximate 
Hamiltonian' — precisely, an effective time-independent interaction picture Hamil- 
tonian such that the associated evolution operator can be explicitly determined - 
as it is done for applying the RWA? 

2. Moreover, does the RWA really allow to obtain a correct first order approximation 
(with respect to the perturbative parameter A) of the evolution operator? 

3. Finally, is it possible to cast in a unique theoretical framework the perturbative 
treatment of time- independent and time-dependent Hamiltonians in such a way that, 
for instance, one can study both the time-dependent Hamiltonian Hq + Hi(£) and the 
time-independent Hamiltonian Hq + H$ using essentially the same approach (and, 
hopefully, finding comparable results)? 

As it will be seen later on, the answers to these questions (in the respective order) are the 
following: 

1. No, it is not necessary. One can apply a time-independent perturbative approach 
that allows to obtain approximate expressions of the evolution operator in the re- 
markable form of a one-parameter group of unitary transformations (see sects. 0] 
andEJ). 

2. No, the RWA does not provide, already at the first perturbative order, the correct 
expression of the evolution operator associated with the Hamiltonian H; 4 it allows 
only to reproduce the qualitative behavior of the correct first order expression (see 
sect. IHJ). 

3. Yes. As it will be shown in the subsequent sections, one can develop a suitable per- 
turbative approach that allows to treat on the same footing both time-independent 
and time-dependent Hamiltonians (see sect. EJ). For instance, one can apply this 
approach to the time-dependent Schrodinger picture Hamiltonian Hq + flj(i) (or 
Hq + H^(t)) and to the time-independent interaction picture Hamiltonian Hq + H$ 
obtaining the same result, i.e. the same perturbative expansion of the (Schrodinger 
picture) evolution operator (see sect.EJ). 

4 It should be clear that this statement concerns the evolution operator itself and does not exclude the 
possibility that the behavior of certain experimentally observable quantities, with specific initial conditions 
of the system, may be rather well predicted using the RWA; see the discussion in sect. |7| 
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3 Basic assumptions and strategy 



Let us consider a time-dependent perturbed Hamiltonian H(X;t), namely a selfadjoint 
linear operator of the form 

H(X;t) = H (t) + flo(A;t) , (14) 

where Ho(t) is a selfadjoint (and, in general, time-dependent) operator — the 'unperturbed 
component' — and H (X;t) is a time-dependent perturbation; precisely, we will assume 
that A i — > H <> (X;t) is (for the perturbative parameter A in a certain neighborhood of zero 
and for any t) a real analytic, selfajoint operator-valued function, with H (O;t) = 0. A 
real analytic function can be extended to a domain in the complex plane. Keeping this 
fact in mind, we will specify that a given property holds for A real. For instance, the 
analytic function A i— > H (X; t) will take values in the selfadjoint operators for A real only. 
Let U(X;t,to) be the evolution operator associated with H(X;t), with initial time to: 

iU(X;t,t ) = H(X;t)U(X;t,t ), U(X; t , t ) = Id , (15) 



where the dot denotes the time derivative and we have set H = 1. Then, we have that 



t/(A;Mo) = U (t,t )T(X;t,to), (16) 

where Uo{t,to) and T(X;t,to) are respectively the evolution operator associated with the 
unperturbed component Ho(t) (evolution operator which, if the unperturbed Hamiltonian 
is time-independent, Ho(t) = Ho, is obviously given by e~ lHo ^~ to ^) and the evolution 
operator associated with the interaction picture Hamiltonian 

H(X; t, t ) := U (t , t) H^X; t) U (t, t ) . (17) 

Let us notice explicitly that, since H(0;t,to) = 0, we have: 



T(0;Mo)=Id. (18) 

We will suppose that the unperturbed evolution Uo(t, to) is explicitly known. Then the 
problem is to determine perturbative expressions of T(X;t,to). To this aim, the central 
point of the paper is the assume for T(X;t,to) the following general decomposition: 

r(A;Mo) = exp(-iZ(A;Mo)) ex V (-i£ C{X;t,t )dtj exp (i Z(X; t , t )) , (19) 

where (X;t) i— > Z(X;t,to), (A;t) i— ► C(X;t,to) are operator-valued functions which depend 
analytically on the perturbative parameter A; in agreement with condition (]18[). we set: 



Z(0;t,t ) = 0, C(0;t,t ) = 0, Vt . (20) 

We stress that the presence of the term exp(i Z(X; to, to)) in formula (|19f) ensures that 
T(X;to,to) = Id, allowing the possibility that Z(X;t, to) ^ for t = to- 
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It will be seen that decomposition (|19|) has a wide range of solutions and that a 
possible choice for fixing a certain class of solutions is given by imposing the condition 
C(X;t,to) = C(A), i.e. assuming that the function (X;t) i— > C(X;t, to) does not depend 
on time. This decomposition includes, as particular cases, two decompositions of the 
evolution operator that have been considered in the literature: 

• the decomposition that is obtained setting Z(X; t, to) = 0, Vi, in formula ()19|) . decom- 
position which is at the root of the Magnus expansion of the evolution operator [S]; 

• the classical Floquet decomposition that holds in the case where the interaction pic- 
ture Hamiltonian depends periodically on time (say with period T) — decomposition 
which is obtained setting C(X;t,to) = C(X), Z(X;to,to) = 0, and assuming that 
(A, t) i — > Z(X; t, to) is periodic with respect to time with period T — and that is at 
the root of the Floquet- Magnus expansion of the evolution operator jllj . 

From this point onwards, for notational convenience, we will fix to = 0. Then, decom- 
position ()19|) can be rewritten as 

T(X;t) =exp(-iZ(X;t)) exp^-i^ C(A;t)dtj exp(iZ(A)), (21) 

where: T(A; t) = T(X;t,0), Z(X;t) = Z(X;t,0), Z(X) = Z(X;0), C(X;t) = C(X;t,0). Let 
us now proceed to obtain perturbative expansions of C(X;t) and Z(X;t). To this aim, if 
we require the interaction picture evolution operator to satisfy the Schrodinger equation, 
we get: 

H{X;t)T(X;t) = if(X;t) 

= e -iZ{X,t) j^sz^t) z(X;t)e- isZ ^)ds e -/oC(A;t)d V Z(A) + 

+ e -iZ(X;t)j^ e -isJ^C(X;t) dt c{x . t) e isf*C(X;t) dt\ dg e -i/*C(A;t) dyZ(A) ^ ^ 

where we have used the remarkable formula (see, for instance, ref. 

^ e F = e F j 1 (e~ sF F e sF ) ds = jf * (e sF F e~ sF ) ds e F , (23) 

which extends to an operator-valued function t >— > F(t) the formula for the derivative of 
the exponential of an ordinary function. Next, let us apply to each member of eq. (1221) 
the operator e lZ ( X '^ on the left and the operator e^ tZ ^e l ^o c ( X '^ dt on the right: 

Ad exp(iZ(A;t)) H(X; i) = J o (Ad cxp(isZ(A;t)) Z (A; t) + Ad cxp( _ is/otc(A;t) dt) C(A; t)) ds , (24) 

where we recall that, given linear operators X, Y, with X invertible, Ad^y := 
Then, since X is of the form e , we can use the well known relation 

oo ^ 

Ad exp(x) y = exp(ad x )^ = ^^ad^y, (25) 

fc=0 
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with ad^ denoting the k-th power (a,d x = Id) of the superoperator adx defined by 
adx Y := [X, Y\. Eventually, applying formula (|25j) to eq. 11241) and performing the in- 
tegrals, we obtain: 

00 -k 00 -k 00 / ' \ k 

£ h a 4 ( A; t) m t)=j2 t^ttv ^)+zZ tStv ad ^(A; t) dt ■ 

k=0 ' k=0 y '' k=0 V ; ' 

(26) 

This equation will be the starting point for the determination of the operator-valued 
functions (A,i) 1— > Z(X;t) and (X,t) ^ C(X;t) at each perturbative order in A, task that 
will be pursued systematically in the next sections. 

4 The time-independent case 

We will first consider the most important special case: the case where the Hamiltonian (|14l) 
does not depend on time. We have at least three good reasons to single out this case and 
to study it at the beginning: 

• to show that in this case — 'the simplest one' — our perturbative decomposition is 
far from being trivial; 

• to highlight the link between our approach and standard perturbation theory for 
linear operators, a link that on the other hand is completely missing for Dyson and 
Magnus expansions; 

• to show how this special case can be regarded as a natural starting point for extending 
the approach to the general time-dependent case ('induction versus deduction'). 

As it will be clear soon, it will be now convenient to set 

Z{\-t):=U (t)Z(\;t)U G (t)\ Z{X) = Z{X; 0) = Z(X) , (27) 

and re-express eq. (|26|) in terms of the transformed operator Z(X;t). To this aim, let us 
first notice that 

Z(X- 1) = Ad w (i(A; t) - i ad z(x . t) # (*)) • (28) 
Besides, given linear operators X, X and Y, with X invertible, one can show easily that 
ad k AdxX Ad x Y = Ad x ad k x Y, k = 0,1,2,... . (29) 

Then, since Z(X;t) = Adjy^t Z(X;t), H{X;t) = Ad^^t H (X;t) and relation 1)28(1 holds, 
using formula (J2H|) . from eq. (|2fi|) we obtain: 

00 -k 00 -k 

Ad t/oWt£^ ad W)#o(A;i) = Ad C7oW t(^-^— ad| (A;4) Z(A;t) 

k=0 ' k=0 ^ 

00 .£> oo / 

" E^ ad W) W)) +E7^^ c(A:t)dt C(A;t). 

k=l ' k=0 V ; ' 
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Next, applying the superoperator Ad[/ (t) to each member of this equation and rearranging 
the terms, we get 

^ a d k z{X;t) (H (t) + H o (X-,t))+H o (X-,t) = Ad^ o(i) ^ 7 ^-ad / fc ic(A;t)dt C(A; i ) 
fc=l * k=0 ^ 

+ E7fcTTv ad l ( A;t)^(A;t). (30) 

k=o y >' 

This equation is, in general, harder to solve than eq. H26|) : but, in the time-independent 
case, we have that H (t) = Hq, ff (A;t) = H (X), and it is natural to assume: 

C(X;t) =C(X) , Z(X;t) = Z(X;0) = Z(X;0) = Z(X) . (31) 

Then, eq. ()30|) can be recast in a much simpler form: 

00 -k 

E % u ^z(X } (Ho + H (X)) + H (X) = e~^ C(X) e iH °\ (32) 

k=l 

Now, observe that the first member of this equation does not depend on time, hence the 
function t \—* e~ lHot C(X) e lHot must be constant. It follows that, if we want eq. ()32j) to be 
consistent, we have to assume also that [C(A), Hq] = 0, i.e. that C(A) is a constant of the 
motion for the unperturbed evolution generated by Hq. Eventually, we find: 

00 -k 

E h a 4(A)(#o + H oW) + ^o(A) = C(A) . (33) 

k=l 

At this point, we are ready to obtain perturbative expansions of the operators C(A) and 
Z(X) (hence, of the interaction picture evolution operator T(A;i)). We will suppose that 
the unperturbed Hamiltonian Hq has a pure point spectrum, while the case where this 
hypothesis is not satisfied is a particular case of the general treatment developed in sect.|51 
We will denote by E\, E 2 , ■ ■ ■ the (possibly degenerate) eigenvalues of Hq and by Pi , P2 , • • • 
the associated eigenprojectors. Since the functions A 1— > H <> (X), X 1— > C(A) and A 1— > Z(X) 
are analytic and Hq(0) = C(0) = Z(0) = 0, we can write: 

00 00 00 

^(A)=^A n iI n , C(A) = ^A"C n , Z{X) = Y,X l Zn- (34) 

n=l n=l n=l 

Now, in order to determine the operators {C n } ne ^ and {Z n } ne ^, let us substitute the 
power expansions ()34j) in eq. ()33|): in correspondence to the various orders in the pertur- 
bative parameter A, we get the following set of conditions: 

d - 1 [Z x , Hq] - H x = , [Ci, H ] = (35) 
C 2 - i [Z 2 ,Hq] + \[Zi, [ZuHq]} - 1 [Z x , H 1 ]-H 2 = , [C 2 , H ] = (36) 
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where we have taken into account the additional constraint [C(X),Hq] = 0. This infinite 
set of equations can be solved recursively and the solution — as it should be expected (we 
will clarify this point soon) — is not unique. The first equation, together with the first 
constraint, determines Z\ up to an operator commuting with Hq and C\ uniquely. Indeed, 
since 

[d , H ] = => C 1 = Y J PmC 1 P m , and \Z X , H ] = ^ (E l - E 3 ) P 3 Z x P t , (37) 
we conclude that 

= J2 P ™ HlP ™> and Z l = Y,P m Z l P m +iY,(Ei-E 3 Y l P 3 H l P l . (38) 

This last equation admits a minimal solution which is obtained by imposing a further 
condition, namely: P m Z\ P m = 0, m = 1, 2, . . . . 

For n > 1, we will adopt an analogous reasoning. Given an operator X, let us set 

71 Am 

g n (X;Z 1 ,...,Z n ) :=J2 — ad z kl ■ ■ ■ &d Zkm X , (39) 

m=l fciH \-km=n 

with n > 1. Then, for n > 2, we can define the operator function 

n-l 

G n (-ffoi • • • j -ffn; #1, • • • , Z n _i) := Q n ^ m {H m ; Z\, . . . , Z n _ m ) — i[Z n , Hq] + il n . (40) 

m=0 

At this point, one can show that the sequence of equations generated by formula (|33|) is 
given by 

C 1 -i[Z 1 ,H ]=H 1 , [C 1 ,H ]=0 
C n -i[Z n ,H ] = G n (H ,...,H n ;Z 1 ,...,Z n _i) , [C n ,# ] = n>2 (41) 



In order to write the general solution of this sequence of equations, it will be convenient 
to introduce a shorthand notation; given a linear operator X, we set: 

(\X\) Ho :=Y, P mXP m , \)X(\ Hq :=X-(\X\)h = Y. P 3 XP ^ ( 42 ) 

m j^l 

1\) X (\]h ■■= *'E(^ " Ejy'PjXPt . (43) 

Now, assume that the first n equations have been solved. Then, the operator function 
G„, + i(i^o, • • • , Hn+i', Z\, . . . , Z n ) is known explicitly and hence 

C n +i = (\G n+ i(H , . . . ,H n+ x; Zi, . . . , Z n )\) H , (44) 
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[Z n+ i,H ] — i |)G n+ i(Uo, . . . ,H n+1 ; Z lt ..., Z n )(\ Ho . (45) 

Again, this last equation determines Z n+ \ up to an arbitrary operator (|Z n +i|)^ o com- 
muting with Hq; in fact, we have: 

Z n+1 = (\Z n+1 \) Hn + [|)G n+1 (tf , ■ ■ ■ , H n+1 ; Z h . . . , Z n )(\] Ha . (46) 



We stress that, in general, the choice of a particular solution for Z n+ \ will also influence 
the form of C n +2, Z n+ 2, .... Thus, we conclude that the sequence of equations defined 
above admits infinite solutions (even in the case where Hq has a non-degenerate spectrum). 
However, there is a unique minimal solution { M C n , M Z„} n6 fj which fulfills the following 
additional condition: 

<rZ n |) Ho = 0, n = l,2,... . (47) 

To clarify the link of our approach with standard perturbation theory for linear op- 
erators, let us recall a few facts (see P |17j). It is possible to show that, under certain 
technical conditions, there exist positive constants ri,T2, ■ ■ ■ and a simply connected neigh- 
borhood I of zero in C such that, for any A 6 2 and m = 1,2,... , one has that: 

1) the following contour integral on the complex plane 

Pm(X) = ^- < f dz (* - ^Wr 1 ( 48 ) 

- where T m is the anticlockwise oriented circle [0, 2-k] B 9 ^ E m + r m e ld around 
the eigenvalue E m — defines a projection (P m (X) 2 = P m (X)), which is an orthogonal 
projection for A G 1 n R, with P m (0) = P m , and I 3 A h P m {\) is an analytic 
operator-valued function; 

2) the range of the projection P m (X) is an invariant subspace for H (A) (but, if the range 
of P m is not 1-dimensional, in general not an eigenspace), hence 

H(X) P m (X) = P m {X) H(X) P m (X) ; (49) 

3) there exists a (non-unique) analytic family A i— ► W(X) of invertible operators such 
that 

p m = w{xy l p m {x) w{x) , w(o) = id (50) 

— with W(X) unitary for A real — which is solution of a Cauchy problem of the type 
i W'(X) = J(X) W(X), W(0) = Id, where the apex denotes the derivative with respect 
to the perturbative parameter and A i— > J(A) is any analytic family of operators - 
selfadjoint for A real — such that 

Y,Pl(\)J(\)P m (\) = ^P4(A)P m (A) (* P m (A)P4(A)P m (A) = *). (51) 

l^m m 
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In standard (Rayleigh-Schrodinger-Kato) perturbation theory, one can obtain the pertur- 
bative corrections to unperturbed eigenvalues and eigenvectors exploiting formula (|48[) and 
a suitable expansion of the resolvent operator (z — H(X))~ 1 (see, for instance, ref. JHj)- I n 
order to recover our previous results, we can use, instead, properties 2) and 3) (compare 
with [2| IH!)- Indeed, let us define the operator 

if(A) := WiX^HiX) W{\) , (52) 

which, for A real, is unitarily equivalent to H(X). Using relations (|49|) and (|5U|) . we find 

77(A) P m = WiXy'HiX) P m (X) W{X) = W(A)" 1 P m (A) H{X) P m {X) W(X) (53) 

and hence: 77(A) P m = P m H(X) P m , m = 1, 2, .... It follows that [H(X),H ] = and 
then we obtain the following important relation: 

[WiXy^iX) W(X) - H , H ]=0. (54) 

Thus, if we set W(X)" 1 H(X)W(X) - H = (7(A), W(X) = exp(-iZ(A)) , and we apply 
relation (|23|) . we find precisely eq. (J22|)- 

Concluding our treatment of the time-independent case, it is worth stressing that, due 
to conditions (|3*T|) . for the overall evolution operator we have: 

U(X;t) = e ~iHot e -iZ(X;t) e -iC(X)t e iZ(X) 

(* Z(X; t) = e~ iHot Z(X; t) e iHot = Z(X) *) = e~ iZW e~ iHot e^W* e iZ{x) 

(* [C(X),H ] = 0*) = e ~iZW e ~i(H +C(\))t e iZ(\)^ ^ 

or, more explicitly, 

U(X;t) = e- iZ ^ ^expH(£ m + \C\ m {X))t)P m e iZ ^, (56) 

rn 

where we have introduced the reduced rank operators 

[C\ m (X):=C(X)P m = P m C{X)P m , m = l,... . (57) 

Notice that the truncations at each perturbative order of the expression (|56|) retain the 
fundamental property of forming one-parameter groups of unitary transformations. To 
make this statement precise, let us introduce the following notation. Given an analytic 
function A i— ► /(A) = J2^Lo A™ /« an( ^ ^ xe d a perturbative order N, we will set 

N 

/ro(A):=5> n /n; (58) 

n=0 

moreover, given another analytic function A h(X), we will set: 

/(A) « h(X) 44 f m (X) = h m (X). (59) 
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Then, for the evolution operator associated with H (A; t) we have that 
\ N 

U{\; t) exp(-i Z M (A)) exp(-i (#o + C [N] (A)) t) exp(i Z [N] (A)) , (60) 

where [C[ N ](X),Hq] = 0. Thus, the N-th order truncation of our perturbative decompo- 
sition of the evolution operator, i.e. the r.h.s. of relation (|6U[). is indeed a one-parameter 
group of unitary operators. 

We conclude this section observing that one can read the array of eqs. (f53|) 'proceeding 
from the bottom to the top', namely, one may assume the decomposition 

U(X;t) = e~ iZ{x) e -*(*o+c(A))t e iZ(X) ^ [C , (A) ( ffo ] = o , 

as a starting point in the case of a time-independent Hamiltonian and induce from this 
case the general decomposition (|21|). This is actually the path that has led the author to 
find the results presented in the paper, extending in a natural way the results previously 
found in the time- independent case (compare with refs. [21 El El El EJ ) • 

5 The general case 

We will now consider eq. ()26|) in its full generality, equation which can be re-written as 

00 -k / i \ 

E h. ad ^;*) ( #(A; t] - jb + i ^ (A; t} ) = £(A; t] ' (61) 

where 

£(A; t) := C( A; t) + £ T^yy dt C(A; t) . (62) 

k=l 

The operator C(A;i) can be recovered from the operator <£(A;i) by means of an order by 
order procedure. Thus, we can solve eq. (|61|) for (£(A; t) up to a given perturbative order 
and obtain the perturbative expansion of C(X;t) truncated at the same order parallely. 
Indeed, if we substitute in eq. (jf)2"|) the power expansions C(X;t) = J2^Li^" ' C n (t) and 
<£(A;t) = Y^n=i ^ n £n(t) , and we single out the various perturbative orders, we conclude 
that the n-th order, which on the l.h.s. is given simply by A n <£n(t), consists on the r.h.s. of 
X n C n (t) plus a function of C\(t), . . . , C n _i(i) and J" * C\(i) dt, . . . , f C n _i(t) dt. Thus, we 
can achieve an order by order solution. Indeed, one finds out that C(A; t) can be obtained 
from (£(A; t) by the following recursive procedure: 

Ci(t) = d(t), 

C„(t) = R n (c 1 (t),...,C n - 1 (ty,J Ci(t) dt,...,^ C n _i(t) dt) +£ n (t), n>2, 

(63) 



14 



where the operator functions R n (. . .) are defined as follows. Given linear operators 

X,X 1 ,...,X n and Y ll . . . ,Y n , 

let us set 

- (—i) m 

n n {X-Y 1 ,...,Y n ):=-Y J -r^r l £ ady kl ■ ■ ■ a d Ykm X , n>l. (64) 

m=l fciH (~fcm=n 

Then, for n > 2, we can define the operator function 

n-1 

Rn(^i> • • • j X n ^i;Yi, . . . , Yn_i) := lZ n - m (X m ; Y\, . . . , Y n _ m ) . (65) 

m=l 

Let us now investigate the perturbative solutions of eq. (|61|) . Substituting the power 
expansions H(\;t) = En=i^H n (t), C(A;t) = En°=l A " £ n(t), Z(A;i) = £„° =1 A" Z n (i), 
we obtain an infinite set of coupled equations that allows to compute order by order the 
operators {<t n (t)} n ^, {Z n (t)} n& ^. In fact, for n > 1, let us set 

F; Z Xl - . - , Z n ) := £ ^ ad ^ • • • ad *J X ~ ^TT J • ( 66 ) 

m=l fciH \-k m =n 

Then we can define G n (Hi(t), . . . , H n (t); Zi(t), . . . , Z n _i(i); Z\(t), . . . , Z n _i(i)) as 

n-1 

£g n - m (^ m (i),Z m (t);Z 1 (t),...,Z„_ m (t)) + i7 n (t), n>2. (67) 

m=l 

With these notations, one can write the sequence of coupled equations which gives a 
perturbative solution of eq. (|61j) as follows: 

Z x {t) = ^i(t)-£i(A;t), 



Z n {t) = G n i y H 1 (t),...,H n (t);Z 1 (t),...,Z n _ 1 (ty,Z 1 (t),...,Z n ^(t)^ 

- <tn(t) , n > 2 , (68) 



As in the time-independent case, this infinite set of equations can be solved recursively. 
Moreover, by virtue of the recursive procedure (|63|) . one can calculate order by order both 
the operators {C n (i)} ne N and {Z n (t)} n€ fq. Indeed, integrating with respect to time each 
equation in the sequence (|68jl and combining the new sequence of equations so obtained 
with the recursive process (|63|) . we find 

z 1 (t) = y^flittj-dtt)) dt+zx, 

Ci(t) = Ci(t), (69) 
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at the the first order, and 



Z 2 (t) = J ^ 2 (H 1 {t),H 2 (t);^i(t);Zi (t)) -£ 2 (t)) dt + Z 2 

= jf f i ad Zl(t) (^(t) - i Zi(i)) + fla(i) - £ 2 (t)) dt + Z 2 , 

C 2 (i) = R 2 (c 1 (t),J Ci(fidi\ +Ca(t) 

= ^ ad / 'C 1 (i)dt C 'i(*) + £ 2(i), (70) 



In general, at the n-th order, for n > 2, we have: 

Z n (t) = (G n (...,^ n (t);...,Z n _i(t);...,Z n _i(t)) -£„(*)) dt+Z n , 

C n (t) = R n r...,C n _i(t);...,^*C n _i(t) dt\ + <tn(t). (71) 

Here, differently from the time-independent case, at each perturbative order we have a 
couple of equations. The solution of the first couple (|oT?|) is obtained by choosing the ar- 
bitrary operator- valued function t t— > <£i(i) and the arbitrary operator Z±; similarly, the 
solution of the n-th couple of equations, for n > 2, involves the previously computed func- 
tions t i— > Ci(i), . .. , i i— > C n _i(t), i i ^ Z\(t), . .. ,t i— > Z n _i(t) and requires the choice of 
the arbitrary operator- valued function t i— ► (£ n (i) and of the arbitrary operator Z ra . T/us 
choice can be fitted according to computational convenience or adapted to specific require- 
ments. 

While the first point (computational advantages) can be best appreciated by means of 
concrete examples — see sect. HO — the second one (adapting the solutions to some require- 
ment) will be illustrated considering two important issues. One of these — to characterize 
the solutions of our perturbative decomposition given a certain class of interaction picture 
Hamiltonians — will be the subject of last part of the present section. The other one is 
the following. In many applications it is customary to study a physical system described 
by an interaction picture Hamiltonian H(X; t) using some effective Hamiltonian E(X; t) 
which is easier to treat, claiming (usually, on the physical ground) that one can in such 
a way achieve a satisfactory description of the system. Now, it is natural to address the 
question: 

what is the relation, at each perturbative order, between the evolution operator 
associated with the effective Hamiltonian and the true evolution operator? 

In order to establish a precise setting for this question, we will first provide an interpreta- 
tion of decomposition (|21|) which sheds light on its meaning. Notice that this decomposi- 
tion can be rewritten as 

T c (X;t) =T z (\;t)- 1 T(X;t)T z (X;0), (72) 
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where 

T c (X;t) := exp^-i^ C(X;t)dtj and T Z (X; t) := exp(-i Z{\; t)) . (73) 

Formula (|72jl can be regarded as a passage to a further 'generalized interaction picture' 
performed on the Hamiltonian H(X;t). Indeed, let us observe that T z (X\t) satisfies the 
equation 

if z (X;t) = 3 L (X;t)T z (X;t) 

= T Z (X; t) 3 R (A; t) , with 3 R (A; t) = T Z (X; t)' 1 3 L (A; t) T Z (X; t) , (74) 

— where 3i_(A;t), 3r(A;£) have the following explicit form: 

00 / -\k 00 -k 

^ = E TTTTV ad l ( A,) 3(A; , 3 R (A; i) = £ tt^Tv ad ^(^) ^ (A; t} (75) 
fc=0 ^ fc=o ^ 

- or, equivalently, i ( J^-T^ 1 ) (A; i) = 3r(A; i) T Z (X; t) -1 . From this relation and eq. I|72|l. 
one finds that 

i7b(A;i) = (r z (A;t)- 1 #(A;t)T z (A;t)-3 R (A;t))r c (A;t). (76) 

Then, by ({75)) and l|76|). one can conclude that eq. (|6*T|) expresses precisely the fact that 
£(A; t) is the transformed Hamiltonian obtained by switching to this 'new interaction 
picture'; namely: 

£(A; t) = T Z (X; t)' 1 (ff(A; t) - 3l(A; i)) T Z (X; t) = T Z (X; t)' 1 H(X; t) T Z (X; t) - 3 R (A; t) . 
It follows that 

/ 00 rt \ 

" 1 (77) 



T C (A; t) = exp Li X n £ C„(i) dtj 



is nothing but the Magnus expansion of the evolution operator associated with the new 
interaction picture Hamiltonian £(A;t). 5 

Hence, coming back to our initial question, it is now clear that, setting £(A;t) = E(X;t), 
decomposition (|21|) provides precisely an order by order comparison between the true 
evolution operator T(X;t) and the effective one Tc(X;t) (notice that only the arbitrary 
constants {Z n } n ^ are still to be fixed in order to determine the decomposition of T(A; t)). 
A remarkable example of such a comparison will be given in sect. El 

Let us now turn to the issue of characterizing a certain class of realizations of our 
perturbative expansion. To classify the whole range of possible solutions and their specific 
properties would be obviously a problem far beyond the scope of the present paper, in 
its full generality. Our aim is to study a wide but coherent class of solutions that is 
particularly relevant for applications. 

J Again, we stress that we use here the term 'interaction picture' in a generalized sense. 
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We will first focus on the important class of solutions which is determined by the condition: 
<£i(i) = £i(0) = <£%,..., Cn(t) = £n(0) = €n, ... Vt. This condition is equivalent to the 
following: 

C 1 (t) = C 1 (0) = C 1 ,...,C n (t) = C n (0) = C n ,... Vt. (78) 

Moreover, if this condition holds, we have: C\ = (£1, . . . ,C n = £n> • • • • Then the solution of 
the first equation — namely Z\({C\, Zi};t) = f H\(t) dt—t C\+Z\ — is fixed by the choice 
of the 'arbitrary constants' C\ and Z\. Inductively, the solution of the n-th equation can 
be achieved by substituting the previously obtained solutions 1 1— > Zi({Ci, Zi}, t), . . . ,t 1— > 
Z n -i({Ck, Zk}^z[; t) — that are fixed by the choice of the arbitrary constants C±, . . . , C n -\ 
and Z±,. . . , Z n _\ — of the first n — 1 equations in the formula 

, Z n -i({Ck, Zk}^zl;t); ... , Z n _i({Cfc, Z k }^~l;tj)di 
n > 2 , (79) 



Z n {{Ck, Zk}f. = i, t) 



G„( 

t C n + Z n , 



which involves the n-th order arbitrary constants C n and Z n . 

We will now consider a specific possible choice of the arbitrary constants {C n , Z n } n ^ 
up to a certain perturbative order N G N. Suppose that the following limits exist (the 
existence and the meaning of these limits will be investigated presently): 



°Ci := lim - [ Hi(t) dt 



OC r 



T 



\jo 



> Zl := lim - / (H^-^d) dt dt 



= ^Br(I^ (t)dt ) dt+ rS' 

; (so) 

°°C N := lim f ± jTg n (...;...;..., ZN-i({°°C fc) "Zfc}^ 1 ; t)) dt) , 

°°Z N := Urn 1 jTQf 6n (...;...;..., ZN-i({°°C fc , ^fc 1 ; t)) dt) dt + \ r °°Cn) . 

Then, one can set {C n = °°C n , Z n = °°Z n }^ =1 . This particular choice has a remarkable 
property. Indeed, in the time-independent case — Ho(t) = Hq and H (X;t) = H (\) - 
there is a precise relation between the solution associated with the arbitrary constants 
{°°C n , °°Z„}„ g N introduced here and the solutions obtained in sect. EJ As we have done 
in that section, we will assume that the unperturbed Hamiltonian Hq has a pure point 
spectrum. At this point, one can prove that: 

1. for any N G N, the limits (jHOJl exist; 
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2. the minimal solution {^C n , M Z n } ne p : j of the sequence of equations (|4"T|) . i.e. the 
solution obtained imposing condition (|47|) . satisfies 

™C n = °°C n , tx3 Z n = °°Z n , VnGN; (81) 

3. the solution of the sequence of equations ((7$)) determined by the arbitrary constants 
{°°C n , °°Z n } riGN is such that 

Z n ({°°C k , °°Z k } n k zl;t) = e lHot °°Z n e~ iHot , Vn G N . (82) 

Notice that the operators {°°C n , °°Z n } ne N can be calculated by formulae (|8U|) where the 
eigenprojectors of Hq — that are involved in the formulae of sect.0] — do not appear. 

We will now prove the existence of the limits (|8Uj) (and provide a simple interpretation 
of their meaning) in two important cases where the (Schrodinger picture) Hamiltonian H 
of the quantum system is, in general, time-dependent: 

• the case where the functions t 1— > Hi(t), . . . are periodic; 

• the case where the functions 1 1— > Hi(t), . . . are operator- valued trigonometric polyno- 
mials — in the following, for the sake of conciseness, just trigonometric polynomials 
— namely functions of time of the type 

m 

F(t) = Y, A ke i " kt , ^el, (83) 

k=l 

where A\, . . . , A m are operators in the Hilbert space of our quantum system. 

Periodic functions and trigonometric polynomials are examples of almost-periodic (a. p.) 
functions (see |191 12U1 l2*T] ). 6 This class of functions has the relevant property of giving 
rise to a Fourier analysis that generalizes the standard Fourier analysis associated with 
periodic functions (a generalization mainly due to the great mathematician Harald Bohr). 
It is worth mentioning that, not only periodic functions, but also trigonometric polyno- 
mials play a prominent role among a. p. functions, since any a. p. function can be suitably 
approximated by a trigonometric polynomial (consider, for instance, the standard Fourier 
expansion of a periodic function). 
Let us recall a few properties that will be useful soon: 

(a) for any a. p. function F — in particular, for any periodic function or trigonometric 
polynomial — one can define the mean value § F § of F, namely it exists the limit 

{ Fl := lim - / F(t) dt; (84) 
T ^°° r J 

6 The general definition of a a.p. function is rather technical and will not be used in the following, so we 
omit it here and we address the reader who may be interested to the cited references. A typical example of 
a C-valued function which is a.p., but not periodic, is provided by the function t i— > (e 1 * 11 + e 1 " 21 ), where 
Xi, ><2 £ R are such that is irrational. 
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if F is a periodic function, with period T, one can easily show that •{[ F ]}- coincides 
with the ordinary mean value of F as a periodic function, i.e. \_F J = ~ f F(t) dt; 
if F is a periodic function with period T (alternatively, a trigonometric polynomial), 
then F — -{[ F ]} is a zero- mean- valued periodic function with period T (respectively, 
a zero- mean- valued trigonometric polynomial); 

(b) given a a.p. function F, if the primitive 

JF : M 9 f ^ / F(t) dt (85) 

is a a.p. function, then §F]j- = 0; in the case where, in particular, F is a periodic 
function (alternatively, a trigonometric polynomial), the primitive JF is periodic 
(respectively, a trigonometric polynomial) if and only if j[ F ]}- =0; 

(c) for any a.p. function F, one can define the essential primitive of F, 7 namely the 
function 

fF(t):=/(F-{F})(*)= /*F(t)dt-t|F}; (86) 

■/ o 

if F is a periodic function with period T, then its essential primitive fF is periodic too, 
with the same period; moreover, the essential primitive of a trigonometric polynomial 
is still a trigonometric polynomial. 

At this point, we will assume that 

up to a certain order N , the coefficients t \— * Hi (£),..., i i—» H^ (t) of the per- 
turbative expansion of the interaction picture Hamiltonian (A; t) \— * H(X; t) are 
periodic functions sharing a common period T, or, alternatively, trigonometric 
polynomials. 

Then, we can show that the limits (|80[) exist. Indeed, if the function H\ is a periodic 
function with period T (alternatively, a trigonometric polynomial), then, by property (c), 
its essential primitive j-Hi is a periodic function with the same period (respectively, a 
trigonometric polynomial). Notice, now, that 

°°Ci = fHi}, 

°°Zi = -{{ {Hi }} = -{{ f(Hi -{Hi})}, (87) 

and 

Zi({°°Ci, °°Zi}; t) = fHi(t) - {{ m }} ; (88) 

hence, t i— > Zi({°°Ci,°°Zi};t) is a (zero-mean-valued) periodic function with period T 

(respectively, a zero-mean- valued trigonometric polynomial). 

Next, if N > 2, by this fact and our initial assumption, we find that the function 

t ^ G 2 {Hi(t),H 2 (t); Zi({^CirZi};t); ZiirCi^Ziht)) (89) 

This term refers to the idea (typical of signal analysis) that j-F contains the essential spectral infor- 
mation about F; i.e. that F can be reconstructed from j-F up to its constant component {F} ('up to a 
constant offset'). 
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is periodic with period T (respectively, a trigonometric polynomial) and 



'C 2 = {{G 2 {...,H 2 (-);...;Z 1 (rC 1 rZ 1 }; (•)))}}; 



(90) 



moreover, by property (c), we have that the essential primitive of the function (|89|) is 
periodic with period T (respectively, a trigonometric polynomial). As a consequence, the 
limit °°Z2, being (up to a minus sign) the mean value of the essential primitive of the 
function (JHHJ), exists and the function t ^ Z2({°°C k , °°Z k }1 =1 ; t) is periodic with period T 
(respectively, a trigonometric polynomial) again. 

For N > 3 — re-iterating the same argument — one finds out that the limits (jSHJ) exist 
up to the order N. In fact, one finds that, for 2 < n < N, the function 



t '-»• G n ( . . . , H n (t); . . . , Z n _i({°°C k,°°Z k }^ = l; t); . . . , Z n _i({°°C k , °°Z k }^ = l ; t)) (91) 



— together with its essential primitive — is periodic with period T (respectively, a trigono- 
metric polynomial). Thus, we have: 



'C n = { G„( H n ( •);...;..., Z n -i({ °C k ,°°Z k }'^ =1 ; (•)))} > 

>Z n = -{ / G n ( . . . ,H n ( •);...;..., Z n ^i({°°C k , °°Zk}kZi'j ( ' ))) }} : 2 < n < N ; (92) 



Z n ({°°C k , °°Z k } k=1 ; t) — (;fG n ( . . . , H n { •);...;..., Z n -i({°°C k , °°Z k } k= l; ( ■ ))))(*) 

- {fG n {...,H n (-) ] ..r,...,Z n „ 1 ({™C k rz k } k Zb(-)))} (93) 



is periodic with period T (respectively, a trigonometric polynomial). 
In conclusion, we have shown that 

if the coefficients t Hi(t) , . . . of the perturbative expansion of the interac- 
tion picture Hamiltonian are periodic functions sharing a common period T 
(alternatively, trigonometric polynomials) up to a certain order N > 1, then 
the limits ()8(J|) exist, and they are given by formulae (|87|) and ()92() ; moreover, 
for any n < H, the n-th order solution t i— ► Z n ({°°C k , °°Zfc}^ =1 ; t), which is 
given by formula (|88[) or (|93|) . is a periodic function with period T (respectively, 
a trigonometric polynomial); 

in addition, observe that we have the further remarkable property (recall the notation (|59|) 
introduced at the end of sect.QJ): 

if the functions t ^ H\{t) , . . . , t \— » H^(t) are periodic, sharing a common 
period T, in the N-th order approximation of the interaction picture evolution 
operator 



moreover, for any 2 < n < N, the n-th order solution 



T(A;t) exp(-iZ [N] (A;i)) exp(-i C [N] (A) t) exp(i Z [N] (A)) 



(94) 



where 



N 




(95) 



n=l 
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N N 

C [N] (A) = ]T A n °°C n , Z m (A) = Z m (A; 0) = £ A" °°Z n , (96) 

n=l n=l 

i/ie function R9t^ exp(— i Z [N] (A; £)) is periodic with period T. 

Actually, one can easily prove (along the lines previously drawn) that — in the case 
where the coefficients 1 1— * Hi(t) , . . . , 1 1— > Hu(t) are periodic functions sharing a common 
period T (alternatively, trigonometric polynomials), if at each perturbative order n < N 
one fixes the operator Z n arbitrarily, and sets C\ = -jf-Hi^ and, for N > 2, 8 

C n = {{ G n ( . . . , H n ( •);•••;••• , Z n ^({C k , Z k }lz{; (•)))}}, 2 < n < N , 

so that Zi({Ci,Zi};t) = ;fHi(i) + Z 1 and, for N > 2, 

Z n ({Ck, Zk}k =1 ; t) = (fG n (. . . , H n { •);•••, Z n _i({Cfc, Z k }^zl; (•));••• ))(*) 

+ Z n , 2 < n < N , (97) 

- the preceding results remain true; namely: the function t t— > Z n ({C k , Z k Yk=\'i *)j f° r 
any n S {1,...,N}, is periodic with period T (respectively, a trigonometric polynomial). 
Thus, our perturbative expansion of the evolution operator can indeed be regarded as a 
generalization of the Floquet-Magnus expansion, which is recovered — in the case where 
the functions t ^ Hx(t), . . . ,t \— * Hf^(t) are periodic, sharing a common period — by 
setting in particular: Z\ = , . . . , Zn = . 



6 The ion trap Hamiltonian (revisited) 

We will now re-consider, in the light of the theory developed in sects. El IU and |SJ the 
quantum Hamiltonian describing a trapped two-level ion interacting with a monochromatic 
laser field (in the Lamb-Dicke regime). Actually — both for 'pedagogical reasons' and for 
a wider comparison with the literature on this subject — with respect to sect. we will 
consider a slightly more general model. As in sect. [21 the Hilbert space of our model is 
Hf® C 2 where H F is the Fock space, namely a infinite-dimensional Hilbert space endowed 
with an orthonormal basis {\n) : n = 0,1,...}, and with the annihilation and creation 
operators a, associated with this basis: a |0) = , a \n) = y/n \n — 1) , n = 1, 2, . . . . 
The Hamiltonian that we will consider in this section is of the form 

H(\;t)=H + H»(\;t), (98) 

where the unperturbed component Hq, as in sect. [21 is given by 9 

Hq = v Id C 2 + - e Id„ F <g> a z , u, e > (99) 

8 One can check recursively that the function 1 1— ► G n ( . . . , H n (t); Z n -i({Ck, ^fc}^Zj; £)) (hence, 
its essential primitive), for 2 < n < N, is periodic with period T (respectively, a trigonometric polynomial). 

9 In this section, with respect to sect. [5] we will adopt a somewhat more formal notation that highlights 
the tensor product structure of our model. 
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- with n and a z denoting respectively the number operator a) a, n\n) = n\n), and the 
effective spin operator (associated with the internal degrees of freedom of the two-level 
ion) o z |±) = ± |±) , |+) = (1, 0) , |— ) = (0, 1) — while the analytic perturbation H (X; t) 
is now defined by (compare with the interaction term H^(t) defined by (J5J)): 

H (X;t) := Xu(e iat [g{n) + e**{a f + a}))® a_ + e~ iat (g{n) + e~ i4, {a f + a\j\®a\ (100) 

with a± = |±) (=f| , <7 + = a_ (as usual), with denoting the 'deformed oscillator opera- 
tor' in TLy defined by 

a f :=af(n) = f(n + l)a, (101) 

and with the functions g, f : {0} UN-*R (notice that, due to (|101j) . we can set, without 
loss of generality, /(0) = 1) and the phase factor e 1 ^ characterizing the model. Observe 
that, since the functions g, f are assumed to be M-valued, we have: 

g(n? = g(n), f(n? = f(n), a] = f(h) o* = a) f(n + 1) . (102) 

Obviously, setting 

g(n) = l, f(n + l)=ri, Vn£{0}UN, = -i , 

we recover the model, with interaction term H^(t), considered in sect. [2 while if we set 
instead e = u, a = and g(n) = 0, f(n + 1) = Vn + T, Vn G {0} U N, = 1, we 
obtain a generalization with 'counter-rotating terms' (i.e. \v cij (g> <x_ and Xva^ <S> o+) of 
the 'Jaynes-Cummings model with intensity-dependent coupling' studied, for instance, in 
ref. 

Another natural choice of the functions g, f is provided by a more accurate approximation 
of the interaction term Hi(t) of the ion trap Hamiltonian, approximation which is adopted 
by some authors (see, for instance, refs. [231 EH EH] ) • Indeed, it turns out that we have: 

H t (t) = Xv(e iat D(-ir])a- + e~ iat D(i V ) a + ) 
= Xv{e iat e-^ l2 (^ Q { m h) 

+ J2(- i v) m (®™fr,n)a m + J m <!> rn (ri-,n))^ ® a_ + h.c.Y (103) 

m=l ' 

where the operator functions &o(rj; h), $ m (r?; n), m = 1, 2, . . . , are defined by 

00 I ■ \2l 00 I ■ \2l 

z=o ^ '•' z=o ^ '-' 

oo /. -.21 

* n (mn) ==E i| ||; ro)| Ni. (105) 

with: 

[n]o = 1 , Hz = n (n - 1)+ • • • (n - I + 1)+ , / > 1 , 
(n — m)+ = for n < m , (n — m)+ = n — m for n > m . (106) 
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Notice that the functions {0} UNBhh &k(i7> n )> fc = 0, 1, . . . , are M-valued. 
In order to derive this result, it is sufficient to observe that 



D(±irj) := exp ^rbi rj (a + 



g-^ 2 / 2 e ±i7?at e ±iva 




1=0 ^ + 1 1 1 1 ' 

from which formulae (|1U4() and (|1U5|) are readily obtained. Notice also that, recalling the 
expression of the generalized Laguerre polynomial 

o;)=e(;+t)^. (108) 



we have: 



<Mr?;n) = L™( V 2 ), m,n = 0,1, 2,... . (109) 

m + n)\ 



Then, in the Lamb-Dicke regime (ij <C 1), we can set (recall that /(0) = 1): 

g{n) = e' v2/2 ^ (r];n) , /(ra + 1) = 77 e^ 2 $i(r/; ra) , Vra G {0} U N , e^ = -i. 

We are not interested here in performing a complete analysis of the dynamics generated 
by the Hamiltonian H(X;t), analysis which would involve a careful study of the different 
regimes associated with specific ranges of the parameters v, e and a. What is of interest 
to us is to fix these parameters in such a way that most of the relevant features of our 
perturbative expansion can be illustrated. We will assume, then, that the parameters 
u, a, e satisfy the resonance condition 

S = e-a = v, (110) 

which is also the most interesting regime from the point of view of applications. With 
condition (|110() . the interaction picture Hamiltonian has the following expression: 

H(X;t) = exp(iH t) H (X;t) exp(-iH t) 

= Xv(e i{a - e)t g(n) <g> a_ + e i *{e< a - v -'* a f ® a_ + e -<(«-*-K)t a \ tg> a J) + /i.e.) 

= Xv(e- iut g{h) ® a- + {eT i2vt a f <8> <r_ + a ] f ® <r_) + /i.e.) . (Ill) 
Notice that it depends periodically on time, with period T = 2tt/u. 

At this point, we recall that there is an almost ubiquitous empirical rule in quantum optics, 
the rotating wave approximation — aimed to drastically simplify the determination of the 
evolution operator — that says: 
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"in order to compute the evolution operator, skip the rapidly oscillating terms 
in the interaction picture Hamiltonian". 

This recipe, applied to our case, would lead us to consider for the interaction picture 
Hamiltonian the following approximation: 

H(X; t)«Ai/ (e** a] ® ct_ + 0/ <g> a+) =: E(X) . (112) 

Clearly, the effective Hamiltonian E(X) — with f(n) = n , n = 0, 1, . . . , and e*^ = — i - 
coincides with the effective Hamiltonian H^t found in sect. [3 We are going to show that 
the 'approximation' (jll2|) does not even produce the correct first order expression of the 
evolution operator. 

To obtain the first order perturbative expression of the evolution operator associated 
with the Hamiltonian H(X;t), we will first exploit the method outlined in sect.[5J To this 
aim, it will be convenient to introduce the analytic function avxp : C — > C, with 

e z — 1 

avxp(z) = for z ^ , avxp(0) = 1 , (H3) 



and to define, for r > 0, the operator 

T d ■■= ~^J^ H(X;t) dt 

= z^^avxp(— iur) g{n) ® <r_ 



+ e^ (avxp(-i2ur) a f <g> cr_ + cfy <g> cr_) + h.c.j . (114) 

From this expression, we have immediately (recall formulae (|80|). and notice that, in this 
case, H(X;t) = XHi(t)) that the operator °°Ci is given by 

°°Ci = hm T Ci = i/ (e^ aj- ® cx_ + e _i * 0/ ® <r+) . (115) 

Then, we can easily obtain the operator °°Z\\ 

°Z X = lim (^T°°Ci-i f t'Cidt 



T^OO 
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= i g(fi) <g> (<t_ - cj + ) + ^ (e^ <8> <r_ - e -< * a) f ® a + J . (116) 

Next, the expression of the operator- valued function 1 1— > Z\(^°C\,°"Z\\t) is found to be 

^i(°°Ci, °°Zi; t) = + t (*Ci - °°Ci) 

= i «/(n) ® (e^"' a_ - e^' a+) 

+ % - (e-*^-® a f ® a. - e i{ - 2vt ~^ a\ tg> <j + ) ; (117) 

hence, it is periodic on time, with the same period T = 2n/v of the operator- valued func- 
tion t i — > H(X;t). The reader may verify that the second order contributions can also be 
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calculated with a rather modest effort. 

We notice explicitly that formula Q115|) shows that the operator A°°Ci coincides with the 
expression of the effective interaction picture Hamiltonian prescribed by the rotating wave 
approximation, but, due to formulae (|116|) and l|117j) . this prescription does not provide 
a correct first order approximate expression of the evolution operator (we will return to 
this point at the end of the section). Besides, observe that, since °°Zi ^ 0, we are not 
recovering the Floquet-Magnus expansion. 

We stress that we could have obtained this result also by looking for the first order cor- 
rection to the evolution operator associated with the effective Hamiltonian E(X) — as 
explained in sect. [5] — with the only additional requirement that the mean value of the 
function t \— > Zi{°°Ci, °°Zi; t) be zero (condition that fixes the arbitrary constant °°Zi). 

Eventually, the unitary operators generated by the selfadjoint operators (|115|) . (|116|) 
and (|117|) can be explicitly computed. Indeed, using the fact that 

(e ld a f ®a±± e~ %0 o) f <g> a^j = (± (a ] f a f ) ® (a T a±) ± (a f a ] f ) <g> (er± a T )j 

= (±l) m {Hnf m |=F> (T\ + f(n + l) 2m |±) (±|) (118) 

— where we have set, for the sake of notational conciseness, 

f (n) := f(n) ^ for n £ N , f (0) = 1 (119) 

— we find easily: 10 
exp(— iX °°Ci t) = expf — iXu(a\ ® <r_ + a f (g) <r+) t 



= cos (Xvf(n) t) tg> |-) (-| + cos (Xuf(n + 1) t) ® |+) (+| 
_ J ^ sin (X, mt) a] e _^ sin(A,f(n + l)t) _ , 

\ f(n) J f(n + 1) 1 

Next, observe that 

,i2ut _t 



iXZi(°°Ci ) °°Zi;t) \g(h)®(e- ivt e ivt v+) A X ( e ^ a /® a -~ 1 " " " , ~ ) 



(121) 



hence, it will be sufficient to compute the separate exponentials that appear on the r.h.s. 
of eq. p21|) . The computation of the first exponential is straightforward: 



=:V m (X;t). (122) 
Then, using again relation (|118|) . we find: 



f »*—e a f ^a +) _ (|Af(n + 1)) <8> |— ) (— | + cos (|Af(n)) ® |+) (+| 

(2i/t-0 sin(|Af(n + l)) 



^-*)^^at®a + )=:^(A;t). (123) 



In the following we will set, by convention, sin(0)/0 = 1. 
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From formulae (|122[) and (|123|) . we can obtain immediately the expression of exp(zA°°Zi); 
indeed: 

exp(a°°Zi) = exp^AZi^Ci, 00 ^!^)) 

~ g~2 * (Of®o-_- a^(g)cr + ) — Ag(n)®(o-_- cr+) 

= tf 2) (A;0)t^ 1) (A;0)t= t( 2) (-A; 0) Tf^-A; 0) . (124) 

The exponentials above provide a simple explicit form of the first order approximation of 
the evolution operator: 

T(A; t) w exp(-iAZi( 00 Ci, 00 Zi;t)) exp(-iA°°Cit) exp(iA°°Zi) 

4 Vfo(A;t) ^)(A;t) exp(-U°°Cit) lfo(-A;0) Vfo(-A;0). (125) 

It is worth stressing the non-trivial fact that the action of this operator on the standard 
basis {|n) <g> |±) : n = 0, 1, . . .} can be easily computed. 

As already observed in sect.|5J a remarkable feature of the time-dependent Hamiltonian 
H(X; t) is that the associated dynamics can be transformed into the dynamics generated 
by a time-independent Hamiltonian by switching to a suitable interaction picture, i.e. the 
interaction picture determined by the reference Hamiltonian ^ a Id Hp ® a z . Indeed, setting 

Rt :=exp(-^aId nF ® a z tj , (126) 
one can define the time-independent 'rotating frame Hamiltonian' 

0(A) := R\ (h(X\ t)-\a Id Hp ® a^j i? t = O + 0<>(A) , (127) 

where (recall relation (J11U)) ): 

00 := uh (g> Id C 2 + - (e - a) Id Wp ® o z = v\h® Id C 2 + - Id Wp (g) cr 2 J , (128) 



0o(A) := Az^ [g(n) + e^(o / + a\)j <g> <r_ + U(n) + e~^( a/ + a})J a+J . (129) 

We notice explicitly that the further interaction picture Hamiltonian 0(A; t), obtained from 
the time- independent Hamiltonian 0(A) taking as reference Hamiltonian the unperturbed 
component 0o, coincides with the 'old' interaction picture Hamiltonian H(X;t); indeed: 

0(A;i) := e mot O (A) e~ if>ot 

= Ai^e - *"*^) ® cj_ + e i4 >{e- i2vt a f ® cr_ + ® cr_) + /i.e.) 

= iT(A;i). (130) 

Hence, we can refer, without ambiguity, to the operators °°Ci, °°Zi, and Zi(°°Ci,°°Zi;t). 
At this point, the reader may easily verify the following facts: 
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• using the results of sect. |IJ one finds that the operators M Ci and ***Z\ coincide 
respectively with the operators °°Ci and °°Zi as given by formulae Q115|) and (|116|) ; 

• the operator-valued function 1 1— > Zi(°°Ci,°°Zi;t) satisfies the relation 

Zi(°°Ci, °°Zi; t) = e m)t °°Zi e - ** * , (131) 
in agreement with what has been stated in sect.|SJ 

The link between the two descriptions — the one associated with the time-dependent 
Hamiltonian H(X;t) and the one associated with the time-independent Hamiltonian Sj(X) 
- is given by: 

= R t e~ iX °° Zl e-i&o+^C^t e *A°°Zi 

= R t e ~ iS}ot (e if)ot e~ iX °° Zl e~ iSi0 *\ e ~ iXXClt e iX °° Zl 

= e ~ iHot e-^M^Cu^Zyj) e -iA°°Cit e ^ 00 z 1 

~ e - iHot T(X;t) = U(X;t). 

We conclude this section with a final comment. Consider the structure of the first 
order expression of the evolution operator generated by Sj(X); namely: 




W(X) <K(\;t) W(A)t 



It is then evident that — since 5K(A; t) coincides with the evolution obtained by the RWA 
— the true first order evolution operator can be expressed as the result of the application 
of a time-independent unitary transformation (associated with the counter-rotating terms) 
to the RWA evolution. Hence, we can say that the RWA provides at least the first order 
qualitative behavior of the evolution operator associated with Sj(X). On the other hand, 
observe that 

exp(-i0(A)t) W(X) exp(-i£ 1) W(X^ exp(-iA°°C7i t) 

= (W(X) exp(-i$j t)W(X^ exp{if) t)) JH(A;t), 

V V ' 

D(A;t) 

where, as it can be easily checked, the operator D (A; t) gives a non-trivial first order con- 
tribution. Hence, the RWA does not provide, already at the first perturbative order, a 
correct approximation of the evolution operator associated with Jo(A). 
It is also worth observing that the RWA and counter-rotating terms play different roles — 
but with the same 'dignity' — in the first order approximation provided by our perturba- 
tive expansion of the evolution operator. 
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7 Conclusions and just a glance to further applications 



In the present paper, we have introduced a perturbative expansion of the evolution opera- 
tor associated with a, in general time-dependent, quantum Hamiltonian that exploits the 
power of the perturbative approach for a twofold purpose: 

• to obtain unitary approximate expressions of the evolution operator, as in the Mag- 
nus expansion, of which it is a generalization; 

• simultaneously, to achieve computational advantages, in the same spirit of standard 
perturbation theory for linear operators. 

On our opinion, it is a remarkable fact that the time-independent perturbative approach 
(essentially, Rayleigh-Schrodinger-Kato perturbation theory) and the time-dependent per- 
turbative approach (Dyson and Magnus expansions) — that are traditionally regarded as 
completely distinct subjects — can be combined to form such a non-trivial blend. This fea- 
ture of our perturbative expansion makes it possible, for instance, to treat closely related 
models — that may be described, depending on the particular case, by a time-independent 
or by a time-dependent Hamiltonian — using the same general method, so that a direct 
comparison of the results obtained is achievable. 

Remarkable examples of applications of the methods outlined in the present paper are 
the applications to quantum optical systems. In particular, applications to the study of 
laser-driven trapped ions (or ions in optical cavities) with various coupling schemes are 
likely to be very fruitful [2] also in view of the possible implementation, by means of 
such devices, of quantum computers For some of these systems, one can switch to 
a suitable interaction picture to obtain a time-independent effective Hamiltonian, as we 
have seen in the example studied in sects. El and El and apply the results of sect. |3J In 
some other cases, this trick does not work and one can apply the more general results 
of sect. a simple example is provided by the Hamiltonian describing a trapped ion in 
interaction with a bichromatic laser field. In all the mentioned cases, the results presented 
here allow to overcome the severe limitations imposed by the rotating wave approximation 
usually adopted in the literature. 

Although a careful criticism of the RWA is not the subject of this paper, it is worth 
mentioning briefly a few facts. As already observed long time ago by Agarwal |27l I28| . 
who studied spontaneous emission effects in two-level ions, the RWA should not be applied 
directly to the Hamiltonian (hence, for obtaining approximate expressions of the evolution 
operator), but only — possibly and with due caution — in the calculation of some exper- 
imentally observable quantities; if one insists in applying the RWA to the Hamiltonian, 
the agreement with the 'true' behavior of the experimentally observable quantities can 
strongly depend both on the choice of the observable itself and on the initial condition 
of the system. The main point to understand is that the behavior of a single observable 
quantity does not contain, in general, all the 'information' encoded in the evolution oper- 
ator. 

Several authors have investigated the contribution of the 'counter-rotating terms' in var- 
ious quantum optical models, contribution which is completely ignored by the RWA; for 
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instance, perturbative corrections to the energy spectrum |29j . and corrections to the time 
evolution by means of path integral [3U] , perturbative EH E3] and numerical jSU ESI 
techniques have been studied. However, as far as we know, our approach is the first sys- 
tematic attempt at taking into account the corrections to the RWA evolution operator 
from a completely general point of view, which allows also to gain a deeper insight in the 
different roles played by the RWA and counter-rotating terms. 

We conclude this section with a final comment. The method introduced in the present 
paper can be extended in a straightforward way to obtain a perturbative expansion of the 
propagator (or 'evolution family', i.e. the generalization of the unitary evolution operator) 
associated with a perturbed evolution equation in a Banach space (a classical reference 
on, in general temporally inhomogeneous or 'non-autonomous', evolution equations is |36j : 
see also ref. |SZI)- Indeed, our method relies on a pure operator approach. This fact 
implies that it could be applied, in particular, to the study of the dynamics of open 
quantum systems |38l 1391 HOI I41j . namely, of quantum systems subject to a coupling with 
an uncontrollable environment or 'bath', coupling that is responsible for the phenomenon 
known as quantum decoherence. Under certain physical conditions, the dynamics of such a 
system can be described by a propagator whose infinitesimal generator is a 'superoperator' 
acting on the density operators, i.e. on the positive unit trace operators in the Hilbert space 
of the quantum system (that form a convex subset of the Banach space of trace class 
operators). In many important cases (see ref. and the rich bibliography therein), the 
'dissipative component' of the superoperator can be considered as (part of) a perturbation, 
hence, the associated propagator may be expanded using our method. 
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